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Abstract 


Formation and evolution of cirrostratus in response to weak, uniform and constant synoptic 
forcing is simulated using a one-dimensional numerical model with explicit microphysics, in which 
the particle size distribution in each grid box is fully resolved. A series of tests of the model 
response to nucleation modes (homogeneous-freezing-only/heterogeneous nucleation) and 
heterogeneous nucleation parameters are performed. In the case studied here, nucleation is first 
activated in the prescribed moist layer. A continuous cloud-top nucleation zone with a depth 
depending on the vertical humidity gradient and one of the nucleation parameters is developed 
afterward. For the heterogeneous nucleation cases, intermittent nucleation zones in the mid-upper 
portion of the cloud form where the relative humidity is on the rise, because existent ice crystals do 
not uptake excess water vapor efficiently, and ice nuclei (IN) are available. Vertical resolution as 
fine as 1 m is required for realistic simulation of the homogeneous-freezing-only scenario, while 
the model resolution requirement is more relaxed in the cases where heterogeneous nucleation 
dominates. Bulk microphysical and optical properties are evaluated and compared. 

t 

Ice particle number flux divergence, which is due to the vertical gradient of the gravity- 
induced particle sedimentation, is constantly and rapidly changing the local ice number 
concentration, even in the nucleation zone. When the depth of the nucleation zone is shallow, 
particle number concentration decreases rapidly as ice particles grow and sediment away from the 
nucleation zone. When the depth of the nucleation zone is large, a region of high ice number 
concentration can be sustained. The depth of nucleation zone is an important parameter to be 
considered in parametric treatments of ice cloud generation. 
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1 Introduction 


The optical depth of cirrus, one of the controlling factors determining its associated net cloud 
radiative forcing, depends on the cloud ice water path (TWP) and effective particle size [e.g., Foot, 
1988]. Some state-of-the-art general circulation models (GCMs) now predict hydrometeor mixin g 
ratios [e.g., £>e/ Genio et al., 1996; Fowler et al., 1996], However, realistic prediction of cirrus 
optical and microphysical properties requires accurate estimation of the number concentration of 
ice particles generated in the nucleation regime. Using parcel model simulations, Kdrcher and 
Lohmann [2002] developed a parameterization scheme for ice particle number concentration via 
homogeneous freezing nucleation (HF) of aerosol particles, and implemented it into the European 
Center HAMburg (ECHAM) GCM [. Lohmann and Kdrcher , 2002] to ex amin e the aerosol effects 
on the ice cloud and earth-atmosphere radiative budgets. Despite advances in parameterization 
schemes of aerosol effects on ice initiation, our fundamental understanding of the evolution of 
synoptically-forced cirrus still lags. Studies based on parcel models are typically not able to provide 
information about the entire cloud from cloud base to cloud top. Moreover, parcel model studies 
usually assume that the ice particles are lifted with the parcel (no particle fallout or fall-in) and that 
there is no exchange of mass or heat with the environment. Neglect of particle fallout/fall-in is 
questionable for weak forcing conditions, where nucleation may last several to more than 1 0 
minutes. Thus, a model of one-dimension (1-D) or higher [e.g., Jensen et al., 1994a, b; 
Khvorostyanov et al., 2001; Sassen et al., 2002] is needed to adequately estimate cloud bulk 
properties over the entire cloud depth. In our study, a 1-D model with an explicit microphysical 
scheme is used to simulate a column of air lifted by a gentle updraft. 

For synoptically-forced cold cirrus, in which the vertical wind speed is less than 0.1 m s' 1 and 
the temperature is less than -40°C, HF of aqueous solution particles, and heterogeneous nucleation 
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via deposition nucleation and/or immersion and contact freezing of aqueous solution particles can 
be activated. The relative humidity with respect to ice (RH;) required to initiate HF is greater than 
the activation RH; required for heterogeneous nucleation. Therefore, unless an air parcel is devoid 
of ice nuclei (IN), it is expected that, as the air parcel ascends and cools, heterogeneous nucleation 
will be activated first. Depending on the magnitude and duration of forcing, and the number 
concentration of activated IN, HF maybe partially or completely suppressed [Lin et ah, 2002a]. For 
synoptically forced cirrus, heterogeneous nucleation may be the dominant nucleation mode and the 
predicted ice crystal number concentration Ni may be significantly reduced in comparison to a 
cloud formed via the HF process due to the presence of heterogeneous IN [Heymsfield and 
Miloshevich, 1 995; DeMott et al., 1997]. 

From the governing equations to predict the amount of ice particles being generated, the 
particle number concentration is determined by the activation RHi and ice particle formation rate of 
the mode N i , given forcing and temperature. Testing a fair range of these two controlling factors 
provides useful insight into the system because the detailed formulation for heterogeneous 
nucleation modes is not understood well. 

In this study, HF of sulfuric acid haze particles is explicitly formulated by adopting the 
effective freezing temperature scheme [Sassen and Dodd , 1988; Heymsfield and Miloshevich, 

1993]. This scheme predicts ice nucleation in close agreement with water activity based on 
homogeneous nucleation theory [Koop et al., 2000]. For the formulation of heterogeneous 
nucleation (HET), we used the exponential form 

Njfi = N 0 exp[P(S l - S')], ( 1 ) 
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where Mn indicates the number concentration of activated IN; N 0 is a small number equal to 10 3 m' 
3 ; S\ and S* are ice supersaturation ratio and activation ice supersaturation ratio, respectively; and 
P controls the IN activation rate above S’ . Parameters p and S* (hereafter, heterogeneous 
nucleation parameters) may be functions of IN and aerosol species, temperature, and updraft speed. 
This formula is the general form of the parameterization scheme for deposition/condensation 
freezing nucleation of Meyers et al. [1992], hereafter referred to as M92. Parcel simulations of 
immersion freezing in sulfuric acid haze particles by DeMott (results can be found in Lin et al. 
[2002a]) indicate that the exponential form may be used to approximate the exact formulations for 
the heterogeneous and homogeneous freezing, as demonstrated in Figure 1 where results from 
parcel model simulations using exact formulations yield a linear relationship between logarithm N{ 
and RHj. In the parcel model simulations presented here, Ni = Mn before HF is activated. This 
approximation (Eq. 1) is easy to implement in any model; however, it does have a weakness. For 
homogeneous freezing nucleation, nucleation does not cease at the peak RHi. The stochastic and 
time-dependent nature of the formulation of homogeneous freezing nucleation differs from that of 
the activation of cloud condensation nuclei. Thus, this approximation could underestimate the 
predicted ice particle number concentration by homogeneous freezing, if it is not carefully 
implemented. Furthermore, the exponential form treats IN as if there is unlimited number of IN. An 
easy way to improve this is to add an upper bound of Si to the equation, as many have done [e.g., 
Meyers et al., 1992; Khvorostyanov et al., 2001]. 

For our nucleation studies, we simulate a long-lived cirrostratus event that was measured with 
the GKSS Raman lidar [Reichardt et al., 1996] above the Esrange research facility (67.9°N, 

21.1°E) near Kiruna, Sweden, on 16-17 January 1997. In the following, the 1-D model is first 
briefly described in Sec. 2, followed by a summary of the lidar measurement and initial conditions 
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in Sec. 3. The HF scenario is discussed in Sec. 4 with an emphasis on the effects of model 
resolution. The results of HET simulation series are presented in Sec. 5 followed by concluding 
remarks (Sec. 6). 

2 1-D Cirrus Model with Explicit Microphysics 

The cirrus model is the 1 -D time-dependent version of a cirrus model developed to simulate 
tropical cirrus anvil outflow [Lin, 1997]. Ice crystals and H2SO4 aerosol particles are grouped into 
bins according to their ice mass and solute mass, respectively, to resolve the evolution of particle 
size distributions (PSDs). For a given bin k in any grid box, the mean H2O mass (/«*) and number 
concentration (Nk) are both predicted. The bin boundaries (xk) follow Xk+i/xk=c, where c is about 1 .5 
and 2.2, and xj is about 3x 10' 15 and 10' 21 kg for ice crystals and aerosols, respectively. The 
nucleation temperature of ice crystals is also tracked. Any given ice bin is further divided into sub- 
bins according to the nucleation temperature for further analysis as in Reichardt et al. [2004]. 

When two ice crystals collide and form an aggregate, the nucleation temperature of the larger 
particle is assigned to the aggregate. For simplicity, ice crystals in all simulations are assumed to be 
columnar as Bailey and Hallet [2004] concluded that, from diffusion chamber experiments, ice 
crystals are mostly polycrystals with columnar characteristics and columns for T < -40°C and RH; > 
110%. We use a power-law relationship between crystal length and mass [Auer and Veal, 1970; 
Heymsjield, 1972; Mitchell and Arno tt, 1994]. 

The governing equation for particle concentration in the k* size bin, Nk, is as follows, 


dN> 


d\np 


- -y(w-v, t )N t + w^N„ L + d 4r W +% Ur.CS) 




a n > 


dt dz v ' ' dz ' dt lgrow,h dt lagg dt lHF ' dt 


where the first and the second RHS (right hand side) terms, respectively, are the particle flux 
divergence in the vertical direction and a compensating term due to anelastic approximation. The 
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other RHS terms, from left to right, are local tendencies due to particle growth, aggregation, 
homogeneous freezing, and heterogeneous nucleation. W is the vertical air velocity, and V t is the ice 
crystal terminal fall speed. Ice crystals initiated via heterogeneous nucleation^ (HET) always start 
from the first size bin (about 1 pm). HF-produced ice crystals may start from other bins consistent 
with the mass of H 2 O in the freezing aerosol; thus, this setup is useful for freezing of large aerosols 
or activated droplets at warmer temperatures (T>~ 45°C). Ice particle fall speed, a function of 
particle size, is calculated using Bohm [ 1 989] . A second order in time and second order in space 
monotonic upstream-biased transport scheme [Allen et al., 1991] is used to treat the particle 
transport in the vertical direction. 

. Time splitting is applied to treat the imposed rate of uplift and local particle growth 
separately. The time step is 5 sec for the advection terms when the grid spacing Az > 4 m and is 
reduced for a smaller grid spacing, while the time step is always 0.1 sec for particle growth. First, 
the local tendencies of moisture and heat due to vertical advection are calculated as is the advective 
tendency of particle concentration on a bin-by-bin basis accounting for fall speed. Then, the 
particle growth in each grid volume is treated in a Lagrangian fashion where the water vapor v 
mixing ratio and potential temperature are changing with time due to the advection tendencies and 
phase change (depositional growth). The imposed advective tendencies of heat and moisture are 
held constant over the advective time interval and the particle concentration is held fixed at the 
value predicted for the end of the advective time interval. 


A hybrid method-of-moment scheme is adopted to treat the terms 


dt 


1 , dN. k 1 , 

I growth and lagg ™ 


Eq. 2 -[Rokicki and Young, 1978; Ochs and Yao, 1978; Chen and Lamb, 1994], It is assumed that 
the PSD inside a. bin is a piecewise-linear function. Instead of integrating the growth and collection 


7 



equations over the PSD inside a given bin, the mean mass of a given bin is used to estimate the 
mean growth rate and interaction with other bins through collection. The piecewise-linear function 
is then used to estimate mass and number of particles transported (via growth, or decay, and 
aggreation) to other bins. This scheme has been demonstrated to produce fairly good simulations of 
the evolution of PSDs [Chen and Lamb , 1994] and is more flexible than strict method-of-moment 
schemes [e.g., Tzivion etal., 1987]. 


The diffiisional growth rate of ice crystals is calculated using electrostatic analogy and 
prolate approximation for the capacitance of columnar ice crystals [Pruppacher and Klett, 1997] 


dm k 

~dt 


RJ 

D'e, 


4 , 4 
+ - — — ( — — 
K'T R V T 



(3) 


In this study, the deposition coefficient is set to 1.0, because existent data suggest that the 
coefficient decreases with increasing temperature [Table 5.5 of Pruppacher and Klett, 1997; 

Haynes et al., 1992; Heymsfield and Miloshevich, 1998] and the predicted N, is not very sensitive to 
the coefficient when its value is between 0.24 and unity [Lin et al., 2002a; Gierens et al., 2003]. 
Radiation effects on ice crystal growth are ignored. The diffiisional growth of aerosol is also 
explicitly calculated. 


Homogeneous freezing of H 2 SO 4 aerosols is computed using the effective freezing 
temperature approximation [Sassen and Dodd, 1988], (J eff = T + AAT m , where T e n and A 4 are 

effective freezing temperature and melting point depression, respectively. A = 2 .) together with a 
fitted HF formula for pure water [Heymsfield and Miloshevich, 1993]. The aerosol is assumed to be 
lognormally distributed with 7/ a = 200 cm' 3 , pure H 2 SO 4 mode radius of 0.02 pm, and a distribution 
width of 2.3. The number concentration of freezing aerosols is not sensitive to the initial haze PSD 
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when the ascent is weak [Jensen et al., 1994b; Karcher and Lohmann, 2002; Lin et al., 2002b], 
Note that, although we do not use Eq. 1 to calculate homogeneous freezing, we will use the 
associated effective (3 and S * forHF (Table 1) for qualitative discussion. 

Heterogeneous nucleation in the upper troposphere is still a subject of active research. Haag 
et ah [2003] has concluded from the interpretation of aircraft measurements during the 
Interhemispheric Differences in Cirrus Properties From Anthropogenic Emissions (INCA) that 
heterogeneous nucleation as well as homogeneous nucleation plays a role in the formation of cirrus 
in the polluted regions of northern hemisphere. The discrepancy between the laboratory measured 
RHj to activate homogeneous freezing and the highest observed upper-tropospheric RH, implies 
that heterogeneous nucleation might be the dominant nucleation mode when the forcing is weak 
[Heymsfield and Miloshevich, 1995], Table 1 shows the tested values of /3 and S* in this study. 
Although the derived /3 and S* in M92 were for temperatures between -7°C and -20°C and Si 
between 0.02 and 0.25, they are adopted to the temperature and Si ranges considered here and are 
expected to represent the most favorable condition for heterogeneous nucleation. A reasonable 
expectation is that the physically appropriate values of (3 and S* lie between M92 and 
homogeneous nucleation values.. The tested values of (3 are close to those inferred from DeMott's 
parcel study of immersion freezing nucleation; /3 ~ 47 and 60 for the -40°C and -60°C cases, 
respectively (Figure 1). 

Tracking IN by implementing a binned IN category is intuitively straightforward but 
expensive in terms of storage space [Khain et al., 2000]. We adopt this IN tracking approach in this 
study. The number concentration of IN that are activated from S^ to S^k+i, IN*, is calculated 
according to (1). We do not impose an upper bound of Si to equation (1) for any of our M92 or 
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HET simulations. The IN bin boundaries follow 5i,*- 1 —0.01 for <Sj ; * < 0.25; Sijr 1 -0.005 for 

> 0.25. The model treats advection of IN* by the vertical wind explicitly. New.ice particles are 
formed by HET when S{ > S\^+\ and IN* are available, i.e. IN* are not depleted due to prior 
nucleation. Thus, considering a volume of air, once Si > S-^+i, all of the IN in bin Enucleate and IN* 
becomes zero. Larger values of Si (S\ > Sij+ 2 ) must be subsequently achieved for further nucleation 
to occur. Alternative approaches are tested and discussed in Section 5. 

3 Simulation Setup 

The numerical simulations of the long-lived cirrostratus near Kiruna, Sweden, on 16-17 
January 1997, serve two purposes. First, we examined nucleation processes in the context of a well- 
observed synoptically-forced cirrus cloud (this paper). Second, we provided detailed microphysical 
data for calculations of cirrus optical properties, which are presented in the companion paper of 
Reichardt et al. [2004], The lidar observations are shown in Figure 2. The cirrus cloud system is 
associated with a strong northwesterly jet stream (75 m s' 1 ) and surface warm front located over 
central Sweden. The cloud started to develop between 8. 6-9.0 km at around 2010 UTC of 16 
January 1997. As the cloud evolved, many embedded structures were present in the non-uniform, 
and sometimes layered, cloud field. In general, the cloud appeared to be more cellular near the 
cloud top, possibly corresponding to a generation zone, with elongated patterns, resembling 
mesoscale fall streaks [Sassen et al., 1989], that sometimes extended from mid-cloud to cloud base. 

f 

In the work presented here, we do not attempt to explicitly model the small-scale cloud structure 
but focus on the general cloud development (the overall cloud top and base heights, and the 30-min 
averaged cloud extinction-coefficient profiles from three time periods [see Reichardt et al, 2004, 
for details]), which are assumed to respond to the synoptic forcing and used to constrain the values 
of uncertain parameters. 
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Local soundings of atmospheric parameters were not available for that day. Instead, 
temperature and pressure profiles (Figure 3) based on the 0000 UTC radiosonde from Lulea, 
Sweden (65.6°N, 22.1°E) are used model input. Lulea is approximately 265 km downwind of 
Kiruna. The Lulea humidity profile is used as a first guess and, as a tunable parameter, 
subsequently adjusted in sensitivity tests. The other tunable parameter is the synoptic scale rate of 
uplift W, ‘which is assumed to be constant with respect to time and height in the model domain (3.8 
to 10.8 km). A considerable number of simulations were performed to obtain a reasonable humidity 
profile (Figure 3) and the corresponding W (5 cm s' 1 ), that could reproduce the overall 
characteristics of the observations (for details about the optimization process see Reichardt et al. 
[2004]). From simulations of W= 3, 5, 7 cm s' 1 together with several different initial humidity 
profiles; the conclusion can be drawn that the latter affects primarily the development of the cloud 
in the first 3 hours; while the evolution of cloud height, cloud depth and IWP throughout the 7-hr 
simulations depends mostly on W. 

4 Homogeneous freezing (HF) scenario 

For an ascending parcel that is devoid of IN and lifted by a constant speed updraft, 
homogeneous nucleation may continue for 10-30 m of vertical displacement once it is activated 
[Lin' et al, 2002a]. Is it necessary to have a model vertical grid size of 20 m or smaller to simulate 
homogeneous nucleation explicitly? To answer this question, vertical grid spacings of 100, 20, 4, 2, 
and 1 m ire tested. The differences between the simulations are dramatic (Figure 4). The 100-m 
and 20-m simulations exhibit a series of nucleation pulses at cloud top, encircled by the white 
curves in the lower left and center panels of Figure 4, which lead to formation of periodic fall 
streaks. In contrast, only a single major fall streak forms in the 4-m, 2-m (not shown), and 1-m 
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runs. These results demonstrate a strong sensitivity to vertical grid resolution and suggest that very 
fine resolution is required to realistically treat the HF process. 

Sassen and Dodd [1989] and Khvorostyanov et al. [2001] have described similar appearing 
pulsation in their models as cycles of nucleation, crystal growth and decreasing moisture (RHi), 
particle fallout, continued adiabatic cooling and increasing RHi, and reinitiation of nucleation. As in 
Sassen and Dodd [1989], the forcing here is constant and continuous. Conceptually, we can view • 
the air column as a vertical stack of air parcels of small but finite sizes where these parcels are 
lifted by the mean wind W. In the absence of vertical mixing, ice particles may only enter from' the 
overlying parcel and may only exit into the parcel beneath via the fall out processes. Onefc 
nucleation becomes active and ice crystals start to form in parcel j, there is no mechanism for the 
particles to travel to the overlying parcel j + 1 (z(j) < z(j+ 1); z varies with time). As a consequence, 
particles in parcel j cannot affect the RH, evolution in parcel j+ 1 . If parcel j is the uppermost parcel 
containing ice and RHi decreases with height above z(j), then continued forcing (lifting) will cause " 
rising RHj in parcel j+1 and ultimately nucleation in that parcel. At fine resolution, this becomes a 
continuous nucleation process at cloud top. In a Eulerian model with vertical motion forcing, once 
crystals are present in a grid volume, vertical advection will quickly transport some to the overlying 
grid volume where they will be regarded, numerically, as dispersed throughout that volume. There, 
they grow by diffusion if RHi is favorable, reducing RHi (or the rate of RH; increase), and . 
consequently delaying or suppressing the onset of nucleation within that grid volume. We conclude 

that the pulsating nucleation process in the 20-m and 100-m runs is non-physical, an artifact of 

. „ . { - ; ; . .. .. 

inadequately resolving the cloud-top nucleation zone. 

Considering the 4-m, 2-m and 1-m simulations, two distinct regions are evident in the time- 
height displays of IWC (ice water content) and N t (Figure 4). The fall streak region develops in 
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response to the initial nucleation episode between simulation time 1-1.5 hr (Figure 4f) in hi ghl y 
supersaturated air over about a 200-m depth and progressively extends downward as ice crystals 
sediment. The fall streak zone is moderately supersaturated (1 15% > RHj > 100%) after simulation 
time t — 4 hr except for the cloud base sublimating zone. Note that the fall streak contains most of 
the ice mass reflecting the crystal growth that occurs as the particles fall through the supersaturated 
layer. We call the cloudy region above the fall streak region as the cloud top region, in which IWC 
is small, supersaturation is large, and Ni increases with altitude. The cloud top region is formed by 
ice crystals generated in the continuous cloud-top nucleation zone, outlined by the white lines in 
Figure 4f. 

Cloud-top height increases with a speed slightly faster than the mean vertical wind speed, 
because the grid boxes above the moist layer achieve HF nucleation threshold sequentially as the 
entire air column ascends (through the grid) and adiabatically cools. The rate at which cloud height 
increases is determined by W, T, P, nucleation threshold, and particularly, the initial vertical 
gradient of RH; above the moist layer. 

Time series of simulated IWP and the column integrated N\ (N path, hereafter, NiP) are 
compared in Figure 5. It is evident that the 1 00-m run underestimates IWP, and that the nucleation 
pulses in the 100-m and 20-m runs result in significant oscillations of NjP. Furthermore, the IWP of 
the 4-m, 2-m and 1-m runs are quite close but begin to diverge over the last two hours, while the 
N;P of the 4-m and 2-m simulations are significantly smaller than in the 1-m run. An increase in 
grid spacing from 1-m to 4-m results in a slight suppression of the peak RHj (as small as 0.2% 
between the 4-m, 2-m and 1-m runs. Figure 6a), but a large decrease in the particle formation rate 
(Figure 6b). This is a consfequence of the extreme sensitivity of the homogeneous freezing 
mechanism to RHj. The profiles of Ni for the simulations, and the number concentration of freezing 
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haze particle, Wa,fe for the 1-m simulation, are shown in Figure 6c. Note that, = N, for a 
parcel simulation whereas such a relationship does not exist between and N\ in a 1 -D 

simulation. The peak Ni is much smaller than its corresponding M. frr indicating that particle , 
divergence is already very active in the nucleation zone, i.e., the reduced value of N\ in the 1-D 
simulation reflects net loss due to particle fall out. The effects of particle fall speed are very evident 
in the N[ profiles, e.g., values of N t at 40 m below the nucleation zone are nearly half as large as in 
the nucleation zone itself in the 1-m simulation. To further illustrate this process, consider that the 
particle residence time, defined as A z/V t , for 10, 20, and 100 pm crystals is about 6, 1.6 and 0.1. 
min, respectively, in the 1-m simulation. 

The significant differences between the 1-m and the 2-m simulations prompt the question of 
whether the grid spacing has to be reduced still further. To verify the 1-m result, an independent 
parcel simulation using the same microphysical module is done for the parcel with an initial 
altitude of 8.8215 km. If subject to the same lift as in the 1-D simulations, this parcel would reach 
9.5265 km at time 235 min, the altitude of peak RH, in Figure 6. The simulated air parcel reaches 
its peak RHi 160.35% at height 9.525 km, when Ni is about 80 L' 1 . At time 235 min (30 second 
later), N is about 130 L' 1 and RHj becomes 160.22%. Shortly thereafter, nucleation ends, and N{ is 
about 160 L' 1 , which is much higher than the corresponding Ni of the 1-m simulation but is 
reasonably close to 1 50 L' 1 of freezing haze particles at the associated height: ' - ' v 

To further understand the cloud top nucleation zone and the cloud top region, we consider 
trajectories through our 1-D grid. Trajectory A (Table 2) represehts air that goes through the cloud- 
top nucleation zone in the 1-m simulation. A parcel following trajectory A enters the nucleation 
zone from above as the cloud top is ascending faster than the vertical air speed. Shown in Figure 7 
are the time-dependent RH.-lOO, N h and the ice crystal formation rate AA interpolated along the ’ 
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trajectory. In the same figure, the corresponding PSDs are shown for the given time (or height) 
without any interpolation. Along trajectory A, nucleation is active for about 8 min (see inset in 
Figure 7a). Although the peak N, is more than 100 L' 1 , iVj decreases rapidly afterward while the 
PSD narrows and the mean size of the remaining crystals grows in response to diffusional growth 
(Figure 7b). The decrease inJVj manifests the particle flux divergence caused by the outgoing flux 
(downward positive) from the bottom exceeding the influx from the top. Therefore, the assumption 
that ice crystals generated in an air parcel would stay together for a considerable time does not 
apply to synoptically forced cirrostratus. As a result of the particle flux divergence, the water vapor 
uptake (per air volume) is reduced. The change of <Sj with respect to time quantifies this efficiency. 
Following Karcher and Lohmann [2002], 

^ = « i r(S 1 tl) + h + a i (S,+l)]^, (4) 

at at 


where a, = — (— — — ) , and a , = . Consequently, when 

1 T R„C„ R/ q s CRT 


W'-p - I '-d 


- — < — + ^ — 5 ice particles at a given location do not uptake excess water vapor 
dt [a 2 + (S- t +1)] 

dS 

efficiently ( — - > 0 ), and vice versa. The above definition of efficiency differs from the definition 
dt 


of phase relaxation time, (q-q s )/( — — ) [e.g., Khvorostyanov and Sassen, 1998], which does not 

dt 


consider the forcing term, i.e., the first RHS term of Eq. 4. Phase relaxation time could be 
misleading, sometimes. For example, a parcel may have a finite phase relaxation time while the 
supersaturation ratio is actually increasing. To further demonstrate this point, the RHi change rate 
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(dRHi /dt) is calculated and plotted using model output of — , Si, T, andP. As shown in Figure 7a, 

dt 

dRHj/dt becomes positive after t = 4.5 hr and is about 2% per hour at the end of simulation. 

Still considering the 1 -m resolution simulation, how is a layer of large iVj sustained for such a 
long time at the top of the fall streak, although its peak value is decreasing with time? Trajectory B 
represents air that, as it ascends, encounters the fall streak descending from above. Recall that, 
similar to the cloud top nucleation zone, there is a cloud base nucleation zone that persists until it 
encounters the descending fall streak (Figure 4f). Similar to trajectory A, more than 100 L' 1 of ice 
particles are present in trajectory B at the end of nucleation (t = 100 min, Figure 7c). The PSD 
indicates that at this time a small number of 1 00 pm-sized particles have entered the parcel from 
above (Figure 7d). JVi increases dramatically and PSD broadens between 100 and 105 min as the 
lower part of the fall streak descends through the parcel. RH* decreases to 102% at 1 10 min and 
remains close to saturation for 25 min until the fall streak passes. The PSD at time 120 min 
represents the PSD of the upper part of the fall streak. This PSD is narrower, but the mean size is 
larger, about 70 pm. Broad PSDs imply that the histories of existent particles are substantially 
different, i.e., different nucleation time and/or altitude, which is attributed to the 200-ni depth of the 
initial nucleation. In the layer of large Ni, the outgoing particle fliix (dowhward positive) from the 
bottom can be compensated by the influx from the top, in contrast to outgoing particle flux 
exceeding in flux from the top near the cloud top nucleation zone/ The RHj change rate becomes 
positive before the entire fall streak passes. 

The evolution of trajectories A and B clearly shows that, within a few minutes of time, ice 
crystals generated in parcel may exit the parcel and be replaced by incoming ice crystals. The 
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single air parcel study of cirrostratus forced by weak forcing is thus limited to the nucleation stage 
and can not be extended to study the evolution of the cloud bulk properties or mass transport. 

In contrast to the cloud top region, the fall streak features greater IWC and N[, and lower ice 
supersaturation. These drastic differences between the two regions are caused by differences in the 
nucleation zone depth (several hundred meters for cloud initiation resulting in the main fall streak 
versus 10 m for the continuing cloud top nucleation zone) indicating that the depth of the 
nucleation zone is an important parameter in predicting cloud properties. In addition, we find that 
particle flux divergence constantly affects the local N{, even in the nucleation zones. It is also 
evident that the relationship between M in the nucleation zone and N t in the cloud ice mass bearing 
zone (lower fall streak) is not straightforward and is a topic worthy of further exploration. Thus, a 
good estimate of the number of crystals generated in the nucleation regions must be accompanied 
by a good estimate of the aforementioned particle flux divergence effect, which depends on the 
depth of nucleation zone, to parameterize cirrostratus bulk properties in mesoscale or large-scale 
models. 

5 Heterogeneous nucleation (HET) scenario 

We now e xamin e the HET process via the simulations listed in Table 1 . For the ranges of S* 
and P studied here, the time series of IWP and NjP (not shown) indicated that the 20-m setup is 
acceptable (<10% difference) for all except Test S201360 . For consistency, all simulations 
presented here used a grid spacing of 4 m. Homogeneous nucleation is allowed, but does not 
become active because the peak ice supersaturation never achieves sufficient magnitude. 
Aggregation is negligible because of the small M in these simulations. 
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A prominent feature of the clouds in the HET runs is the periodic fall streaks (Figure 8). This 
contrasts with the occurrence of only one dominant fall streak in the 1-m HF simulation shown in 
Figure 4c, f. Layering of the consecutive fall streaks takes place. In addition to the response to the 
initial moist layer, two other nucleation zones emerge: the quasi-continuous cloud-top nucleation . 
zone (small discontinuities due to IN bin discretization), as in the fine resolution HF runs, and an 
intermittent mid-cloud nucleation zone. 

The depth of the cloud-top nucleation zone and number of particles produced Vary among the 
cases (Figure 9). Test M92 (smallest (3 ) has the thickest cloud-top nucleation zone with a depth 
greater than 100 m, the depth is about 25-35 m in other HET tests, compared to 10 m in the HF run 
(largest effective (3 ). This trend indicates that the depth of cloud-top nucleation decreases with /? . 

The peak Si at the cloud top in Test M92 is much greater than the threshold S' . With small (3 
(Table 1), the particle formation rate is small and increases rather slowly as Si increases beyond the 
activation S' . When the particle formation rate is small, depositional growth of the relatively 
small quantity of ice crystals does not uptake much water vapor before they fall out, i.e., the few 
crystals grow rapidly. As Si and the associated particle formation rate become large, the particle 
residence time and number concentration rise, i.e., more particles but smaller sizes. Ultimately, 
depositional growth of the newly activated particles is sufficient to drive dS/dt negative and 
eventually terminate further cloud top nucleation. For the rest of the HET simulations, for a given 
S' within the cloud-top nucleation zone, an increase in (3 is associated with a greater peak particle 
formation rate but smaller peak RHj. Given /i , lowering the nucleation threshold S* also results in 
reduced peak RHi, but leads to an increased particle formation rate and a slight elevation of the 
nucleation zone. The effects of S' and f3 are also reflected in the column-integrated particle 
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formation rates in the cloud-top zone (Figure 10a), where the relative differences are maintained 
throughout the simulations. 

RHj in the HET simulations does not decrease significantly below the cloud-top nucleation 
zone (Figure 9a). This indicates that the particle growth there is insufficient to uptake enough 
water vapor to limit the values of S{. We see from the REf contours in Figure 8 that the upper part 
of the cloud is highly supersaturated. 

Mid-cloud nucleation, which is activated when relative humidity is on the rise and IN are 
available, occurs and acts to limit Si. The nucleation pulse mechanism, briefly discussed in the 
previous section, explains the formation of the mid-cloud nucleation zones. The premise of the 
mid-cloud growth-sedimentation-nucleation cycle is that the cloud-top nucleation zone does not 
produce sufficient number of ice crystals that would sediment into the mid-cloud and limit RHi via 
growth. Below the cloud-top nucleation zone, when particles coming from above cannot effectively 
uptake water vapor, S\ increases in response to continuous forcing ( W ). New ice crystals form when 
IN* are available and their associated S^+i is achieved. The new ice crystals grow and fallout, and 
another cycle begins. Considering a trajectory that passes through the cloud top nucleation zone 
(the first 2.5 hour of the trajectory displayed in Figure 1 1) and then mid-cloud nucleation zones, the 
continuous forcing via W results in S\ achieving not only the corresponding nucleation threshold 
S* (and initial crystal nucleation), but also, later, new successively larger peaks of S{ in the mid- 
cloud region (the trajectory after 2.5 hour in Figure 11) and thus renewed nucleation. In contrast, 
for the HF simulation with high S* and ft , the mid-cloud S\ does not achieve S* and a mid-cloud 
generation zone never forms during the 7-hr simulation. Mid-cloud nucleation acts to broaden the 
PSDs, as illustrated in Figure 12. The PSDs of the HET scenarios are much broader than those in 
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the HF scenario (e.g., Figure 7). The PSDs can be segmented due to the intermittent nature of mid- 
cloud nucleation. •; ... ■. . . i 

We further define cloud-top nucleation as any nucleation taking place in air initially located 
in the layer above 8.2 km before its first peak For example, dRFf/dt is positive initially and 
becomes zero at t = 2.5 hr in the trajectory shown in Figure 1 1 . That is, the air along the trajectory 
is in the cloud-top nucleation zone in the first 2.5 hr. We apply this method to all the trajectories 
above 8.2 km and then establish the lower boundary of the cloud j top nucleation zone as a function 
of time for each simulation. The relative contributions of the cloud-top and mid-cloud nucleation 
for time period t = 1 to 6 hr (i.e., excluding the initial cloud formation event) are plotted in Figure 
1 0. The column-integrated cloud-top contribution decreases with time in response to the decrease 
in the depth of the cloud-top nucleation zone, while its relative contribution, though oscillating, 
also decreases with time. The cumulative cloud-top contribution in the 5 HET cases (t = 2.5 to 6 
hr) is less than 50% (Table 3), reflecting the importance of mid-cloud nucleation when there are 
sufficient IN such that HET is the dominant nucleation process. Ignoring mid-cloud nucleation will 
result in an overestimate of mid-cloud water vapor and an underestimate of cloud ice. 

The IWP, N;P, and optical depth of all the simulations (4 m for all HET cases, and 1 m for the 
HF case) are shown in Figure 13. Nucleation starts later in the HF simulation because more time is 
required for the forcing (W) to raise $ to the required threshold, which is much larger than for HET. 
Thereafter, the HF simulation quickly produces the greater values of NiP until t = 6 hr. The values 
of IWP and optical depth in HF also exceed those in the HET cases, but not till after 3 hours. This 
implies that if the large-scale forcing does not sustain for a long time, the HET scenario might be 
more efficient in producing ice and transporting water downward via particle sedimentation. The 
grouping of HET results (Figure 13) indicates that these bulk properties are more sensitive to S* 
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than ft (for 0.2 < S*'< 0.3 , and 40 < fi < 60) although larger /? does entail a slight increase of JVj 
and optical depth. The properties of the M92 simulation are close to the S20 cases, although 
significant differences exist in the morphology of the cloud (Figure 8). 

Figure 13 also shows that the trend of optical depth follows that of IWP but is further 
enhanced by differences in the bulk mass extinction coefficient, defined as optical depth divided by 
IWP, which depends on the effective particle size of the entire cloud layer. Furthermore, in this 
type of cirrus with high ice supersaturation and low ice particle count, an increase in Np entails 
'more effective water vapor uptake, and thus greater IWP. This suggests that one cannot treat the 
parameterization of IWP and cloud bulk effective radius as two separate issues for synoptically 
forced cirrostratus. 

Finally, sensitivity tests on alternative diagnostic schemes are presented in comparison to the 
IN tracking scheme we have used here. The schemes diagnose activation of IN by checking (1) 

: cLSi/df >0 if $ > 0 (similar idea may be applied to formula based on temperature and height) and (2) 
Max (IVinCS; )- N { , 0) if S{ > 0, respectively. The former (diagnose 1 in Figure 14) treats IN as if 
they are self-replenishing; the latter (diagnose 2 in Figure 14) treated IN as if they move with the 
* ice particles, not with the air. Both diagnostic schemes overestimate iVjP (Figure 14). The former 
diagnostic scheme does a good job for the cloud-top nucleation zone, but overestimates the mid- 
cloud nucleation and IWP. The latter diagnostic scheme overestimates cloud-top nucleation, but 
underestimates the cloud top RHi as well as mid-cloud nucleation. However, the sensitivity is 
smaller compared to the sensitivity to (j and S* (Figure 14 vs. Figure 13). Therefore, IN tracking 
is not a priority to simulate cirrostratus forced by a smooth and uniform ascent. 
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6 Conclusion 


A one-dimensional cirrus model with explicit microphysics is used to study the formation arid 
evolution of a cirrostratus forced by weak vertical motion. A series of tests on the nucleation modes 
(homogeneous-freezing-only versus heterogeneous nucleation) and heterogeneous nucleation 
parameters are performed. The study shows that, when heterogeneous nucleation is the dominant ... 
nucleation process, the complex structure of cirrostratus may result solely from the episodic nature 
nucleation, particle growth and sedimentation in die absence of turbulence and cloud-scale 
dynamical processes. In particular, in order to accurately predict the number of ice.crystals and • 
their size distribution within cirrus cloud layers, we find that the depth of nucleation zone, usually . 
much much smaller than the vertical grid spacing for the upper troposphere in GCMs and cloud 
resolving models (CRMs), should be explicitly represented in parametric treatments of ice cloud 

| ' jif;. 

generation. 

Other important findings are: 

• Three types of nucleation regions exist in the cloud studied. First, the nucleation 
responds to the prescribed moist layer and occurs over a depth defined by the initial 
humidity structure. This would correspond to cloud initiation (formation). Second, in 
the presence of continuous adiabatic forcing, a continuous nucleation region occurs at 
cloud top and elevates with time with a speed greater than the imposed gentle updraft. 
This region is caused by haze particles or ice nuclei reaching freezing/nucleation 
threshold in the RHi-decreasing-with-height layer that bounds the initially prescribed 
moist layer from above. Contribution of the cloud-top nucleation zone has been 
ignored by parcel simulation studies. Last, for the heterogeneous nucleation cases, 
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intermittent nucleation occurs in mid-cloud regions in which ice crystals are formed 
when RHj is on the rise and IN are available. 

The depth of the cloud-top nucleation zone depends on P , given background 
thermodynamic variable profiles. For the HF only scenario (large effective P ), the 
cloud-top nucleation zone is very thin and a model resolution as fine as 1 m is 
required to obtain a good approximation of the generated N{. For the HET scenario, a 
resolution of 20 m may be required. 

Intermittent mid-cloud nucleation reduces water vapor and broadens the PSDs in the 
cloud. A combination of large S* and extremely large ft (the HF scenario) results in 
complete disappearance of the mid-cloud nucleation zone in this study. In contrast, 
small P and, especially, a small S* entails a greater contribution by mid-cloud 
nucleation. The contribution of mid-cloud HET nucleation is more sensitive to S * , 
than to P , for the ranges studied here, and do mi nates the overall particle production 
by the cloud. Systematic differences in PSDs exist between the homogeneous- 
freezing-only (narrower PSDs) and heterogeneous nucleation (broader PSDs) 
scenarios in the studied case, due to the mid-cloud nucleation and greater depth of 
nucleation zones in the heterogeneous nucleation scenario. 

In all simulations, the mid-upper portion of the cloud is quite humid (RH > 120%) 
except within the initial fall streak in the HF scenario. With regard to in-cloud ice 
supersaturation, the cirrus studied here does not represent a rare case. High ice 
supersaturation has been frequently observed by the ARM Raman lidar at the 
Southern Great Plain (SGP) site [Comstock et ah , 2004] and also in 
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subtropcial/tropical anvil and tropical cold thin cirrus [ Gao et al., 2004]. In the low 
ice particle count cirrostratus presented in this study, IWP and excess water vapor 
depend on ice particle path NjP, i.e., more ice crystals enhance the uptake of water 
vapor and thus also the cloud IWP. 


7 Notation 


C k 

D' 

g 
K' 
L 
L s 
m k 
No 

N a 

Ni 


N, 


Njn 


Capacitance for bin k 

Specific heat of air at constant pressure 

Modified diffusivity of water vapor in air 

Saturation vapor pressure over plane ice surface 

Acceleration of gravity 

Modified thermal conductivity 

Ice particle length 

Latent heat of sublimation 

H 2 O mass for bin k 

10 3 m' 3 

Haze particle number concentration 
Ice crystal number concentration , ( 

Ice particle formation rate 
Activated IN number concentration 
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p 

q 


q* 


R a 


R v 

RHi 

£ 

S- 

T 


Teff 


A T 


V t 


w 

A z 

P 

e 

P 


Pressure 

Water vapor mixing ratio 
Saturation water vapor mixing ratio 
Specific gas constant for dry air 
Specific gas constant for water vapor 
Relative humidity with respect to ice 
Ice supersaturation ratio 
IN activation ice supersaturation ratio 
Temperature 

Effective freezing temperature 
Melting point depression 
Terminal fall speed 
Synoptic scale vertical wind 
Vertical grid spacing 
IN parameterization coefficient 
0.622 

Air density 


Subscripts k denote the bin number. 
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10 Figure Captions 

Figure 1 . Ice particle number concentration JV, in an ascending parcel (4 cm s' 1 ) as a function of the 
parcel RHj. Curves Wa004 and Ca004 (both HF and immersion freezing modes are allowed) are 
from parcel simulations using DeMott's model with starting temperatures -40°C and -60°C, 
respectively [see Lin et al., 2002a, for details]. For reference, the result of Ch004 (HF only) is also 
plotted. 

Figure 2. Height versus tim e display of backscatter coefficient generated from consecutive lidar 
profiles of the cirrus cloud measured on 1 6- 1 7 January 1 997 over northern Sweden. For 
experimental details, see Reichardt et al. [2002]. 

Figure 3. Radiosonde data (temperature is indicated by the solid curve and potential temperature is 
denoted by the dash curve) taken at Lulea at 00 UTC, on 17 January 1997. The RHi profile, which 
is denoted by the dotted curve, is selected after several sensitivity tests. 

Figure 4. Effect of vertical grid spacing on the simulated IWC (mg m' 3 ) and N x (liter' 1 ) for the HF 
scenario. From left to right, Az = 100, 20, and 1 m, respectively. The white curves in the upper 
panels indicate RHi (dotted: 60 and 80%, solid: 100, 115, 125, 135%), while the white curves in the 
lower panels outline the nucleation zones. 

Figure 5. Effect of vertical grid spacing on IWP and NiP. 

Figure 6. Effect of vertical grid spacing on the cloud-top nucleation zone of the HF scenario at 
simulation time t- 235 min. Note the thin solid curve in panel (c) is the number concentration of 
freezing haze particles from the 1-m simulation. 

Figure 7. The evolution of trajectories A (panels a and b) and B (panels c and d). RHi-100 (solid), 
jV,. (dotted), AT (dash), and dRHjAfr (dash-dotted) are displayed in panels (a) and (c). Particle size 
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distributions (PSDs) are displayed in panels (b) and (d), where corresponding time is indicated by. 
the corresponding markers in the inset panels in (a) and (c). i •• - , 

Figure 8 . Similar to Figure 4, but effect of heterogeneous nucleation parameters on IWC (mg m’ ) 
and (liter' 1 ). From left to right: M92, S20B60 , and S30660. . ■? , >•„ ; 

Figure 9. Effect of heterogeneous nucleation parameters on cloud-top nucleation zones at 
simulation time / = 3 hr. Note that the lower segments of N- curves (z < 9.3: km) for. cages ; < 

S20B40 and S3 06 60 indicate mid-cloud nucleation. 

Figure 10. Column-integrated particle formation rates for the cloud top nucleation zone (a) and for 
the entire air column (b) after the initial nucleation burst. The relative contribution of cloud-top 
nucleation is plotted in (c). 

Figure 1 1 . The evolution of a parcel trajectory released at zo- 8.1)66 km in Test M92,' Where RHj- 
100 is indicated by the solid curve, M is denoted by thddbtteii curve, and dRHi/dt is indicated by 
the dash-dotted curve. 

• • ' ' . ■ ‘ '■ U l , : d . . v ' i.: . ' •' .'1 

Figure 12. Left panels: Ni and IWC profiles of Test M92 at simulation time t = 5 hr. Right panels: 
PSDs at given heights, indicated by the corresponding markers in the left panels. 

Figure 13. A comparison of the HF and HET bulk microphysical properties (IWP and NiP) and 
optical properties at wavelength 0.55 pm (optical depth and mass extinction coefficient of the 
cloudy layer). Optical depths were calculated with the optical-nncrophysicai cirrus model' of 
Reichardt et al. [2004] assuming ice particles of exclusively columnar hexagonal shape. 

Figure 14. Effect of IN tracking/diagnosing schemes on IWP and NiP. 
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11 Table 


Table 1. Simulation identifiers. M92 adopts coefficients from Meyers et al. [1992], 


Nucleaton Mode 

Run ID 

S* x 100 

fi 


M92 

~5 

~13 


S20640 

20 

40 

HET 

S20640 

30 

40 


S20660 

20 

60 


S30B60 

. 30 

60 

HF 


~59.5 a 

~1 12 a 


a For quick comparison purposes only. These values are derived for the parcel released at z=8.822 
km. See Sec. 4 for details. 


Table 2. Trajectories A and B 


Parcel ID 

^0 

Nucleation Starting 



Time (min) 

Height (km) 

A 

8.663 

126 

9.041 

B 

8.2 

95 

8.485 
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Table 3. Relative contribution of the cloud-top nucleation between simulation time t -2.5 and 6 hr. 


Run ID 

top/all 

M92 

38.7% 

S20660 

34.4% 

S20B40 

40.0% 

S30B60 

48.2% 

S30B40 

45,7% 



RH. (%) 


Figure 1. Ice particle number concentration N\ in an ascending air parcel (4 cm s -1 ) 
as a function of the parcel RHj. Curves Wa004 and Ca004 (both HF and immersion 
freezing modes are allowed) are from parcel simulations using DeMott’s model with starting 
temperatures -40°C and -60°C, respectively [see Lin et al, 2002a, for details]. For reference, 
the result of Ch004 (HF only) is also plotted. 
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Figure 2 Height versus time display of backscatter coefficient generated frbffi consecutive 
lidar profiles of the cirrus cloud measured on 16-17 January 1997 over northern Sweden. 
For experimental details, see Reichardt et al. [2002]. 
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Lulea at 00 UTC on 17 January, 1997 



Figure 3. Radiosonde data (temperature is indicated by the solid curve and potential 
temperature is denoted by the dash curve) taken at Lulea at 00 UTC, on 17 January 1997. 
The RHi profile, which is denoted by the dotted curve, is selected after several sensitivity 
tests. 
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Figure 4. Eff ect of vertical grid spacing on the simulated IWC (mg m ) and N- t (liter ) 
for the HF scenario. From left to right, Az = 100, 20, and 1 m, respectively. The white 
curves in the upper panels indicate RH, (dotted: 60 and 80%, solid: 100, 115, 125, 135%), 
while the white curves in the lower panels outline the nucleation zones. 
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Figure 6. Effect of vertical grid spacing on the cloud-top nucleation; zone, of the HF, scenario ./ 

at simulation time t = 235 min. Note the thin solid curve .in panel (c) is the number 
concentration of freezing haze particles from the 1 m simulation. ...... 







Figure 7. The evolution of trajectories A (panels a and b) and B (panels c and d). RH;- 
100 (solid), (dotted), iVj (dash), and dRHj/dt (dash-dotted) are displayed in panels 
(a) and (c). Particle size distributions (PSDs) are displayed in panels (b) and (d), where 
corresponding time is indicated by the corresponding markers in the inset panels in (a) and 
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Figure 8. .Similar to Figure 4, but effects of het^rpgeneoug, nucleation jjgxametprp on JW,U 
(mg m -3 ) and N\ (liter -1 ). From left to right:, M92, S2Q/360, and S30/360. , , , 
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Figure d. Eff ect of heterogeneous nucleation parameters on cloud-top nucleation zones at 
simulation time t — 3 hr. Note that the lower segments of iVj curves (z j 9.3 km) for cases 
S20/340 and S30/360 indicate mid-cloud nucleation. 
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Figure 10. Column-integrated particle formation rates for the cloud top nucleation zone (a) 
and for the entire air column (b) after the initial nucleation burst. The relative contribution 


of cloud-top nucleation is plotted in (c). 
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Figure 11. The evolution of a parcel trajectory released at zq — 8.666 km in Test M92, 
where RH;- 100 is indicated by the solid curve, TV; is denoted by the dotted curve, and 
dRHj/dt is indicated by the dash-dotted curve. 
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Figure 12. Left panels: N\ and IWC profiles of Test M92 at simulation time t = 5 hr. Right 
panels: PSDs at given heights, indicated by the corresponding markers in the left panels. 
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Figure 13. A comparison of the HF and HET bulk microphysical properties (IWP and NiP) 
and optical properties at wavelength 0.55 /rm (optical depth and mass extinction coefficient 
of the cloudy layer). Optical depths were calculated with the optical-microphysical cirrus 
model of Reichardt et al. [2004] assuming ice particles of exclusively columnar hexagonal 
shape. 
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Figure 14 . Effect of IN tracking/ diagnosing schemes on IWP and NjP 
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Large-scale cirrus cloud systems are an important component of the climate system due 
to their radiative effects. Proper account of such systems in climate system models is 
therefore essential. These models are using increasingly sophisticated treatments 
including representations of ice crystal size distributions and nucleation processes. This 
work provides guidance on the manner in which such parameterizations should be 
formulated in order to accurately account for the effects of nucleation on ice crystal 
number densities and particle sizes, and ice water path. Formation and evolution of 
cirrostratus in response to weak, uniform and constant synoptic forcing is here simulated 
using a one-dimensional numerical model with explicit microphysics, in which the 
particle size distribution in each grid box is fully resolved. The response to nucleation 
modes (homogeneous-ffeezing-only versus heterogeneous nucleation) and heterogeneous 
nucleation parameters are performed. A continuous cloud-top nucleation zone with a 
depth depending on the vertical humidity gradient develops. For the heterogeneous 
nucleation cases, intermittent nucleation zones in the mid-upper portion of the cloud also 
form where the relative humidity is on the rise, because existent ice crystals do not uptake 
excess water vapor efficiently, and ice nuclei (IN) are available. Vertical resolution as 
fine as 1 m is required for realistic simulation of the homogeneous-freezing-only 
scenario, while the model resolution requirement is more relaxed in the cases where 
heterogeneous nucleation dominates. Bulk microphysical and optical properties are 
evaluated and compared. Ice particle number flux divergence, which is due to the vertical 
gradient of the gravity-induced particle sedimentation, is constantly and rapidly changing 
the local ice number concentration, even in the nucleation zone. When the depth of the 
nucleation zone is shallow, particle number concentration decreases rapidly as ice 
particles grow and sediment away from the nucleation zone. When the depth of the 
nucleation zone is large, a region of high ice number concentration can be sustained. The 
depth of nucleation zone is an important parameter to be considered in parametric 
treatments of ice cloud generation. 


